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We investigate the ground-state property and the excitation gap of a one-dimensional spin-orbital 
model with isotropic spin and anisotropic orbital exchange interactions, which represents the strong- 
coupling limit of a two-orbital Hubbard model including the Hund's rule coupling (J) at quarter 
filling, by using a density-matrix renormalization group method. At J—0, spin and orbital correla- 
tions coincide with each other with a peak at q=n/2, corresponding to the SU(4) singlet state. On 
the other hand, spin and orbital states change in a different way due to the Hund's rule coupling. 
With increasing J, the peak position of orbital correlation changes to q=n, while that of spin corre- 
lation remains at q—n/2. In addition, orbital dimer correlation becomes robust in comparison with 
spin dimer correlation, suggesting that quantum orbital fluctuation is enhanced by the Hund's rule 
coupling. Accordingly, a relatively large orbital gap opens in comparison with a spin gap, and the 
system is described by an effective spin system on the background of the orbital dimer state. 

PACS numbers: 75.45.+j, 75.10.Pq, 75.40.Mg 



I. INTRODUCTION 

It has been recognized that orbital degree of freedom 
plays a crucial role in the emergence of exotic magnetism 
in strongly correlated electron systems. The interplay of 
spin and orbital degrees of freedom triggers off a variety 
of orbital ordering that stabilizes a certain spin structure 
on the orbital-ordered background, leading to complex 
ordered phases, as frequently observed in transition metal 
oxides and /-electron compounds^ 

In addition to the appearance of diverse spin-orbital or- 
dered states, the combined quantum effects of spin and 
orbital degrees of freedom yields a possibility of novel 
quantum critical and quantum disordered states^ In this 
context, a high symmetric SU(4) exchange model, which 
describes a spin- 1/2 system coupled with a pseudospin- 
1/2 for two-fold orbital degeneracy, has been one of the 
intriguing subjects from a theoretical viewpoint li^&L^ 
In such a high symmetric case, spin, orbital, and com- 
bined spin-orbital degrees of freedom play an identical 
role£ In particular, the one-dimensional model is Bethe 
ansatz solvable, 4 and the quantum critical behavior of an 
SU(4) singlet state has been revealed by analytical and 
numerical investigations^^ In fact, correlation func- 
tions show power-law decay with a critical exponent 3/2, 
and the elementary excitation is gapless. 

The SU(4) spin-orbital exchange model represents the 
strong-coupling limit of a two-orbital Hubbard model at 
quarter filling, in which electrons hop between the same 
types of orbitals with equal amplitude and the Hund's 
rule coupling is ignored. Indeed, such a simplified model 
is an appropriate starting point to gain deep insight into 
the complex quantum fluctuations of spin and orbital de- 
grees of freedom. However, in actual materials, the high- 
est SU(4) symmetry is likely to be broken down, thus it is 
quite important to clarify the effects of possible symme- 
try breaking. For instance, intensive studies have been 
devoted to a one-dimensional SU(2) sp i n x SU(2) or bitai ex- 
change model and the ground-state phase diagram has 



been revealed ,M24Li2J£i&i£ Around the SU(4) point, a 
quantum critical phase extends in one side next to the 
SU(4) point , 13 i 14 i 15 while in the other side there occurs a 
gapped phase, where spin- and orbital-singlet dimers are 
formed in an alternating pattern^ 

In a more realistic situation, however, the Hund's rule 
coupling causes anisotropic exchange interactions in the 
orbital part, and the SU(4) symmetry is broken down to 
SU(2) spin xU(l)orbitai^> i£ ^> i2 To clarify the effects of 
such symmetry breaking, an SU(4) Hubbard model per- 
turbed by the Hund's rule coupling has been analyzed by 
means of renormalization-group and bosonization meth- 
ods*^ It has been proposed that the excitation gap opens 
for an arbitrarily small Hund's rule coupling, while the 
opening excitation gap is exponentially small as a func- 
tion of the Hund's rule coupling, corresponding to a gen- 
eralized type of Kosterlitz-Thouless transition.— 

On the other hand, it has been common understanding 
that the Hund's rule coupling in multi-orbital systems 
plays an important role in itinerant ferromagnetism i 21 : 22 
In fact, it has been shown that the Hund's rule coupling- 
induces a ground-state transition from a paramagnetic 
(PM) state to a ferromagnetic (FM) state in the strong- 
coupling region ; 23 : 24 ' 25 but there, the discussion has been 
focused on the appearance of ferromagnetism itself and 
the property of spin and orbital states in the PM phase 
has not been clear yet. Thus, we believe that it is an 
intriguing issue to clarify how the SU(4) singlet state is 
deformed to a spin- and orbital-singlet state due to the 
Hund's rule coupling. 

In this paper, we investigate the ground-state property 
and the excitation gap of a one-dimensional spin-orbital 
model with the SU(2) sp i n xU(l) or bitai symmetry, by using 
numerical techniques. When the Hund's rule coupling is 
zero, the system has the SU(4) symmetry, and spin and 
orbital correlations have a peak at q—Ti/2. On the other 
hand, the Hund's rule coupling induces FM and antifcrro- 
orbital (AFO) interactions, leading to the characteristic 
change of the peak position of orbital correlation to q=n, 
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while that of spin correlation remains at q=ir/2 in the 
PM phase. Moreover, we observe the stabilization of an 
orbital dimer state with a significant orbital excitation 
gap. Taking into account the robust orbital dimerization, 
the low-energy physics is described by an effective spin 
system on the background of the orbital dimer state. 

The organization of this paper is as follows. In Sec. II, 
we describe the model Hamiltonian and the numerical 
method. Starting from a two-orbital Hubbard model, a 
spin-orbital exchange model in the strong-coupling limit 
is introduced. In Sec. Ill, we show our numerical results 
of physical quantities such as correlation functions and 
excitation gaps. We discuss distinctive changes of spin 
and orbital states from the SU(4) singlet state due to the 
Hund's rule coupling. Finally, we summarize the paper 
in Sec. IV. 



w -f + e w 4- -i- ir 4- - 4^P 

-f- + « -f ^ + ^ 4^ -« 

FIG. 1: Spin-orbital configuration in (a) one-electron state 
and (b) two-electron state. 

TABLE I: Eigenenergy and eigenstate of two-electron state 
in the ascending order of eigenenergy, where |t<t) denotes the 
spin a and orbital r state. 

eigenenergy eigenstate 

cTW \a W T>, l« I), 7f (l« T)l/3 l)+\<* DI/3 T» 

U'+J 7y(M)|/U>-M>|0T» 

U-J' 73(M>MH/3t}|/U» 

U+J' -75(M>M>+|/3t>l/U» 



II. MODEL AND METHOD 

Let us consider doubly degenerate orbitals on a one- 
dimensional chain with N sites including one electron per 
site (quarter filling). The two-orbital Hubbard model is 
described by 

H = t ^2 {4 Ta d l+lTa + h.c.) + U ^2 p ir1 p lTl 

+ J> 4rArl d "'l d "'h C 1 ) 

where di Ta (d\ Ta ) is an annihilation (creation) operator 
for an electron with spin <r(=|, J,) in orbital t(=oi, (3) at 
site i, PiTcr=d ira .di T(T , and t is the hopping amplitude. 
Here, we assume that electrons hop between the same 
types of orbitals with equal amplitude. Note that, in gen- 
eral, the hopping amplitudes of multi-orbital systems are 
evaluated from the overlap integral between orbitalsi 26 ' 27 
The present simple form of the hopping amplitudes rep- 
resents the system with, for instance, {p Xl p y ) orbitals or 
(d yz , d zx ) orbitals on a chain along the z axis. 

In the interaction terms, U, U', J, and J' denote intra- 
orbital Coulomb, inter-orbital Coulomb, inter-orbital ex- 
change (Hund's rule coupling), and inter-orbital pair hop- 
ping interactions, respectively. These Coulomb integrals 
are all positive. Note that the spin diagonal part of the 
Hund's rule coupling causes an attractive interaction, so 
that the inter-orbital Coulomb interaction is effectively 
reduced to U'—J. To ensure that the total inter-orbital 
Coulomb interaction is repulsive, we consider the region 
of J<U'. We also note that the relation U=U'+J+J' 
holds due to the rotational invariance in the local or- 
bital space, and J=J' due to the reality of the orbital 
wavef unction^ When the Hund's rule coupling is zero, 
i.e., U=U' and J=J'=0, the system possesses the SU(4) 
symmetry^ 



In order to consider an effective model in the strong- 
coupling limit, as usual, we start from the ground state 
of the atomic limit i=0, and treat the electron hopping 
as a perturbation! 28 ' 29 In the ground state of the atomic 
limit, each site is occupied by one electron, leading to 
four- fold degenerate states, as shown in Fig. 1(a). Then, 
the electron hopping causes virtual two-electron states, as 
shown in Fig. 1(b). There appear four eigenenergy states, 
as listed in Table I. Among them, the lowest energy state 
is spin triplet, and has energy U' — J. Thus, the Hund's 
rule coupling stabilizes spin polarized states. The other 
three states are spin singlet. Note that the second- and 
third-lowest energy states are degenerate because of the 
relation U =U'+J+J'. Note also that all these four levels 
are degenerate when U=U' and J=J'=0. 

Now we apply the degenerate perturbation theory by 
considering the second-order processes of the electron 
hopping, and obtain an effective spin-orbital exchange 
model, such as 

i 

+B (s, ■ S, - (T, • Tj - 2T?T? + ~) 
+C (Si ■ S j - ^ (t, ■ Tj 2T?T* + 1) 
+D (s t ■ Sj (t« • T S - 2TVTV + ±) (2) 

where 

S, = (1/2) 4ra°W<W (3) 

is the S=l/2 spin operator and 

T, = (1/2) 4 ra cT TT ,d tT , a (4) 
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indicates the T=l/2 pseudospin operator representing 
two orbitals, with the Pauli matrix er. Concerning the ex- 
change interactions A^D in Eq. ((2]), each of them arises 
from the corresponding eigenenergy of the virtual two- 
electron state, and they are expressed as 



A = it 2 / ({/' — J), (5) 

B = tt 2 /(U' + J), (6) 

C = 4t 2 /(U — J') = B, (7) 

D = At 2 /{U + J 1 ) = At 2 /(U' + 3J). (8) 



Note that the relation A>B>D>0 holds in the realistic 
parameter region of 0<J<U'. Hereafter, we take t 2 /U' 
as the energy unit. Normalized by t 2 /U', the exchange 
interactions are parametrized by J/U'. Since we focus 
on the region of 0< J<U', we set U'=l for simplicity and 
consider the variation of J within 0< J<1. 

Here, let us discuss the symmetry of the spin-orbital 
model. When U=U' and J=J'=0, the spin-orbital model 
is reduced to 

»V = W (*-% + i)(T.-T, + i), (9) 

which has the SU(4) symmetry!^ On the other hand, in 
the spin-orbital model of Eq. @, we can clearly see that 
exchange interactions of the spin part are isotropic, while 
those of the orbital part show anisotropic forms according 
to the energy level of the virtual two-electron state. Note 
that the second and third terms have the same energy 
denominator due to the relation U=U'+J+J' , so that 
exchange interactions in the orbital part become isotropic 
in the zx plane. In fact, Ttot = Si commutes with H , 
and becomes a good quantum number. Concerning spin 
degree of freedom, total spin 5t t and <S? t = Si Sf are 
good quantum numbers. Thus, the SU(4) symmetry is 
broken down to SU(2) sp i n xU(l) rbitai due to the Hund's 
rule coupling. 

In order to discuss the change of the spin-orbital state 
from the SU(4) symmetric case, it is useful to rewrite the 
Hamiltonian by the sum of the SU(4) symmetric part and 
the other parts, such as 

H id = K{Si ■ S s + 1/4) (Ti • Tj + 1/4) 

+J s S t ■ Sj + J Z T X {T*T* + T[-T 3 Z ) + JJfflTj 



+j y ST s, ■ s;/ ;« /;/. (10) 

where the exchange interactions are given by 

K = A + D, (11) 

J s - (—A + B)/2, (12) 

Jff = (A-D)/2, (13) 

4 = (A-B)/2, (14) 

J| T = 2{B-D). (15) 



We can easily see that Js is negative, while J^ x and Jj, 
are positive, indicating that FM and AFO interactions 




J J 

FIG. 2: Exchange interactions in the spin-orbital model (I10|) : 

(a) J s /K, (b) J} x /K, J V T IK, (c) J y ST /K, and (d) K, as a 
function of J. 



are induced by the Hund's rule coupling. We note the 
relation of 

jv + J| T /4 = Jf, x , (16) 

among anisotropic orbital exchanges. 

In Figs. 2(a)-(d), exchange interactions are plotted as a 
function of J. At J=0, exchange interactions other than 
K are equal to zero, and the system possesses the SU(4) 
symmetry. On the other hand, the Hund's rule coupling 
induces FM and AFO interactions, which should cause 
the deviation of the spin-orbital state from the SU(4) 
symmetric case, as will be discussed in detail in the next 
section. 

We analyze the spin-orbital model (flTJ|) . or equiva- 
lently ([2]), by a density-matrix renormalization group 
(DMRG) methodi 30 ! 31 The finite-system algorithm is em- 
ployed with open boundary conditions. We keep up to 
800 states in the renormalization step and the truncation 
error is kept around 5xl0~ 6 or smaller. We also use a 
Lanczos method for the analysis of a four-site periodic 
chain. 



III. RESULTS 

A. Four-site chain 

To grasp the fundamental property of the spin-orbital 
configuration in the ground state, first we consider a four- 
site periodic chain, which is a minimal model to form a 
unique SU(4) singlet ground state at J=04 The SU(4) 
singlet state for N=4 is composed of different four spin- 
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FIG. 3: (a)-(c) Three classes of relevant spin-orbital configu- 
rations in the ground state of the four-site system, (d) Weight 
of each class of Figs. 3(a)-(c) in the ground state as a function 
of J. 

orbital configurations at four sites, expressed as 

isu(4)> = (1/V24) £ dLAaApAiM' 

where |0) is the vacuum state and the summation is taken 
over all permutations in terms of the site index. Note that 
|SU(4)) consists of 24 states with the same weight. For 
finite J, however, the ground state is not represented by 
|SU(4)) itself. In fact, the ground state is changed to a 
spin- and orbital-singlet state, which is well described by 
the superposition of the 24 states included in |SU(4)), but 
the 24 states are split into three classes according to the 
weight in the ground state ) 18 ! 19 as shown in Figs. 3(a)-(c). 
Note that each class has eight equivalent states due to the 
translation and the reversal of orbitals. In Fig. 3(d), we 
show the J dependence of the weight of each class m in 
the ground state, defined by 

^ m = El^G>| 2 , ( 18 ) 

where -0G is the ground state and cf>i denotes the basis. At 
J=0, each class contributes to the ground state with the 
equal weight due to the SU(4) singlet ground state. On 
the other hand, with increasing </, w a gradually increases, 
while Wb does not change so much and w c decreases. Note 
that the total weight of the three classes remains almost 
unity even when J is increased, but it abruptly decreases 
at around J=0.4, since the ground state turns from the 
singlet state into a FM state. 

Here, it is worth noting that each of the three classes 
of Figs. 3(a)-(c) is identified by a distinct peak position 
in spin and orbital structure factors, denoted by (g s , q ): 
(tt/2, 7t) for the class (a), (tt/2, ir/2) for the class (b), and 
(it, tt/2) for the class (c). In the class (a), a parallel spin 
arrangement is stabilized by the FM interaction due to 
the Hund's rule coupling, while the pairs of parallel spins 
point opposite directions to form a spin-singlet ground 
state for small J. Concerning the orbital state, an AFO 
configuration appears so as to avoid the energy loss due 
to the double occupancy in the virtual hopping process. 




distance r 

FIG. 4: (a) Spin and orbital correlation functions at J=0. 
(b) Spin and (c) orbital correlation functions at J—Q.2. 

On the other hand, the class (b) also exhibits a parallel 
spin arrangement, but there occurs a partly ferro-orbital 
configuration, leading to relatively small Wb comparing 
with u> a . The class (c) is more unfavorable due to a fully 
antifcrromagnetic configuration. 

B. DMRG results 

From the discussion within the small four-site chain in 
the previous subsection, we can intuitively understand 
the effect of the Hund's rule coupling to cause the char- 
acteristic change of the spin-orbital state, but it is highly 
required to clarify what types of correlations develop in 
the thermodynamic limit. For this purpose, here we move 
on to our DMRG results for longer chains. In Fig. 4, we 
show spin and orbital correlation functions, defined by 

C„(r) - (1/N r ) ( S ? S j)' (!9) 

|i— j\=r 

C (r) = (1/N r ) £ <7?T;}, (20) 

\i-j\=r 

where N r is the number of site pairs (i,j) with r=\i— j\, 
and (• • •) denotes the expectation value. Note that we 
average over all pairs of sites separated by distance r in 
order to reduce the effect of open boundaries. At J=0, 
as shown in Fig. 4(a), C s and C exhibit exactly the 
same behavior with a four-site periodicity which is in 
agreement with the previous investigations for the SU(4) 
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FIG. 5: (a) Spin and (b) orbital structure factors at sev- 
eral values of J. (c) Spin and (d) orbital structure factors at 
q=Tv/2 and q—TT as a function of J. 

exchange model.— i For finite J, however, C s and C are 
not equivalent any longer due to the breakdown of the 
SU(4) symmetry. In Fig. 4(b), it is found that C s keeps 
the four-site periodicity, while the amplitude of the os- 
cillation becomes large. On the other hand, as shown in 
Fig. 4(c), we observe a drastic change in the structure of 
C . Indeed, the four-site periodicity almost disappears, 
and the correlation decays much faster for long distances 
in comparison with that at J=0. 

In order to see clearly the changes in the spin and 
orbital structures, it is useful to measure spin and orbital 
structure factors, defined by 

%) = T,( s 3 s ® eiqU ~ k) / N > w 

j,k 

T z (q) = J2( T j T k)e lqU ~ k) /N, (22) 

j,k 

T y (q) = J2( T j T k) eiq(J ~ k) l N - ( 23 ) 

j,k 

In Figs. 5(a) and 5(b), we show S(q) and T z (q) at several 
values of J, respectively. At J=0, S(q) and T z (q) agree 
with each other and have a peak at q=ixj1 due to the 
clear four-site periodicity. Even when J is increased, the 
peak position of S{q) remains at q—n/2, since S(n/2) is 
enhanced and S(tt) is suppressed. The J dependence of 
S(tt/2) and S(tt) is plotted in Fig. 5(c). On the other 
hand, the peak position oiT z {q) changes from q=n/2 to 
q=ir, since T z (n/2) is suppressed and T z (ir) is enhanced, 



FIG. 6: Log-log plot of (a) spin and (b) orbital correlation 
functions at several values of J. The results are shown for 
distances of four multiples. 



in sharp contrast to the case of S{q). In fact, we find that 
T z (ir) becomes larger than T z (w/2) at around J=0.05, 
as shown in Fig. 5(d). It is interesting to note that these 
characteristic changes of the spin and orbital structure 
factors are naturally understood from the discussion for 
the four-site chain in the previous subsection. Namely, 
the spin-orbital configuration in Fig. 3(a) is stabilized 
by the Hund's rule coupling at least for short distances, 
leading to the enhancement of S(w/2) and T z (tt) in the 
whole system. 

We note here that T v (q) shows qualitatively the same 
behavior with T z (q) , but there appears a slight difference 
between them. In Fig. 5(d), T z (q) and T*(g) at q=n/2 
and q=n are plotted as a function of J. With increasing 
J, T v (ir/2) is suppressed and T y (w) is enhanced in the 
same way as T z (q), while T z (n) is much enhanced than 
T v (tt). This anisotropic behavior is naturally explained 
by the anisotropy in the orbital exchanges (J-p, Jf., Jg^) 
within a mean-field discussion. For the J| T term of the 
spin-orbital model (fTUj) . replace Si ■ Sj by (S, • Sj). In 
the fully polarized FM state, (Si • Sj) takes the maximum 
value 1/4 at all bonds, and the anisotropy disappears, 
since the relation JlL + Jg T /A = J-f? holds, as shown in 
Eq. (fTB]) . However, in the spin-singlet state, (Sj • Sj) 
should take a smaller value than 1/4, indicating that the 
orbital correlation in the zx plane becomes larger than 
that along the y direction. 

In Figs. 6(a) and 6(b), we show C s and C in the log-log 
scale to observe the decaying behavior. At J=0, C s and 
C exhibit power-law decay with the critical exponent 
3/2, consistent with the previous analytical works^ and 
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FIG. 7: (a) Spin and (b) orbital dimer correlation functions, 
measured from the center of the chain, at several values of J. 
(c) Spin-orbital configuration in the dimer state. Dotted and 
solid oval enclosures represent spin- and orbital-singlet pairs, 
respectively. 



numerical investigations. 5,7 With increasing J, as shown 
in Fig. 6(a), C s gradually increases in accordance with 
the enhancement of S{tt/2), while the slope seems to be 
gentle. In particular, we observe that C s decays as 1/r 
at J=0.2, as denoted by the dotted line in Fig. 6(a). 
At J=0.3, C s comes to decay exponentially, indicating 
the behavior of a quantum disordered spin state. On 
the other hand, as shown in Fig. 6(b), C turns to show 
an exponential decay for smaller J in comparison with 
C s , suggesting that a quantum disordered orbital state is 
stabilized by the Hund's rule coupling. 

However, it is important to remark that, according 
to the generalized Kosterlitz-Thouless transition ; 1 ^ 20 the 
correlation length would diverge as £^exp(const./ J) for 
small J . Unfortunately, from numerical calculations, it 
is quite difficult to distinguish whether the correlation 
length is finite but very large or certainly infinite. In 
this sense, our results indicate that the spin correlation 
develops at least within the present system size, and the 
variation of the power-law decay could be a crossover 
phenomenon from the SU(4) critical state to a quantum 
disordered spin state due to the finite size effect. 

Concerning the order parameter at finite J, it has been 
proposed that there occurs a dimer order with the alter- 
nating pattern of spin- and orbital-singlet pairs P which 
is similar to the dimer state of the SU(2) spin xSU(2) or bitai 
exchange models In order to clarify how spin and orbital 
dimerizations develop due to the Hund's rule coupling, we 
investigate spin and orbital dimer correlation functions, 
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FIG. 8: (a) Spin and (b) orbital correlation functions between 
nearest-neighbor sites at J=0.2. 



defined by 



A>(*)i) — S i+1 (Sj Sj +1 Sj +1 Sj +2 )) 7 (24) 
D„(i,j) = (T?T* +1 (T?T? +1 -T? +1 T* +2 )). (25) 



In Figs. 7(a) and 7(b), D s and D a are shown, respectively, 
for several values of J. Note that, taking into account the 
alternation of the spin and orbital dimerizations, we shift 
the starting point of D s and D Q by one, i.e., we set i—48 
for D s and i—47 for D Q with 96 sites. At J=0, we find no 
indication of dimer long-range order, since both D s and 
D rapidly decay as the distance becomes large. With 
increasing J, D s and D show a tendency to grow in the 
system, but we find a clear difference between them. It is 
observed that D s remains small even when J is increased, 
while D Q shows a significant development, indicating that 
the orbital dimerization is stabilized by the Hund's rule 
coupling. 

In Fig. 7(c), we schematically depict the spin-orbital 
configuration in the dimer state. Here, we note again that 
the Hund's rule coupling induces FM and AFO interac- 
tions. Indeed, with increasing J, the ground state turns 
to be a fully polarized FM state and the orbital state is 
described by a T—l/2 AFO critical state. However, the 
ground state remains spin singlet for small J. It is con- 
sidered that orbital-singlet dimers are formed to stabilize 
a parallel spin configuration in each orbital-singlet dimer, 
while the correlations between the dimers are suppressed. 
Then, there occurs a staggered configuration of parallel 
spin pairs, leading to the spin-singlet ground state. Note 
that such a spin configuration is already seen in the four- 
site problem, as shown in Fig. 3(a). 

We note that the dimer order is associated with a spon- 
taneously broken translational symmetry, and the ground 
state has two-fold degeneracy with a finite excitation gap. 
The broken symmetry is discrete rather than continuous 
rotational symmetry, so that a true long-range order can 
exist even in one dimension. On the other hand, in the 
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FIG. 9: System size dependence of (a) spin and (b) orbital 
gaps at several values of J. Dotted curve denotes the analyt- 
ical result at J=0, given by Eq. (130[l , 



Indeed, it is observed that A s and A Q at J=0 are well 
reproduced by Eq. (j50)) , as denoted by the dotted curve. 
With increasing J, we find that A s becomes small for 
each system size up to 2V=192, as shown in Fig. 9(a). 
Even though the ground state is the gapped dimer state, 
the spin gap does not show a remarkable development, 
since the formation of the spin dimer state is very weak. 
This result is consistent with the analytical work stating 
that the spin gap in the thermodynamic limit is expo- 
nentially small as A s ^exp(— const./ J) 

On the other hand, as shown in Fig. 9(b), we observe 
a quite different behavior of A Q . At the beginning of the 
increase of J (J<0.1), A Q keeps almost the same finite 
size spectrum. For larger J, A D shows a drastic increase 
for each system size, and A Q seems to converge to a finite 
value in the thermodynamic limit. We mention that such 
a development of the orbital gap is accompanied by the 
stabilization of the orbital dimer state due to the Hund's 
rule coupling. Thus, the low-energy physics is described 
by spin excitations below the orbital gap. We conclude 
that the Hund's rule coupling leads to the characteristic 
difference between spin and orbital states, and the mag- 
netic property is represented by an effective spin system 
on the background of the robust orbital dimer state. 



present calculations, one of the two possible dimerization 
patterns is favored due to open boundaries. In fact, we 
can see characteristic spatial patterns in spin and orbital 
correlation functions between nearest-neighbor sites, 

CTW = WSUi), (26) 
C™(i) = (7TI7 +1 >. (27) 

As shown in Fig. 8, we notice that C" n and C" n exhibit an 
oscillation in an anti-phase manner. In particular, truly 
FM/AFO bonds appear around the edges. As shown 
in Fig. 8(b), orbital-singlet dimers are robustly formed 
at the edges and the whole system is efficiently covered 
with the orbital-singlet dimers, so as to lower the energy 
making use of quantum orbital fluctuation. Then, weakly 
formed spin-singlet dimers are arranged alternately to the 
orbital-singlet dimers, as shown in Fig. 7(c). 

Now we turn our attention to the low-lying excitation. 
To clarify how the distinct behavior of spin and orbital 
states is reflected in the spin and orbital excitations, we 
investigate spin and orbital gaps, 

A s = £(1,0) -£(0,0), (28) 
A = £(0,1) -£(0,0), (29) 

where E(S* ot , Tf ot ) denotes the lowest energy in the sub- 
space with S* ot an d ^tof I n Figs. 9(a) and 9(b), we show 
the system size dependence of A s and A Q , respectively. 
We should note that, at the SU(4) point J—0, the elemen- 
tary excitation is gapless, while the finite size spectrum 
includes a logarithmic correction^ 

£ k-ihi)- (30) 



IV. SUMMARY 

In this paper, we have studied the effect of the Hund's 
rule coupling on the SU(4) spin-orbital exchange model. 
The Hund's rule coupling introduces isotropic FM and 
anisotropic AFO interactions as well as the spin-orbital 
coupled term, leading to the symmetry breakdown from 
SU(4) to SU(2) spin xU(l) orb itai- When the Hund's rule 
coupling is zero, spin and orbital correlations coincide 
with each other and have a peak at q=n/2 due to the 
SU(4) singlet ground state. With increasing the Hund's 
rule coupling, the peak position of orbital correlation is 
found to change to q=ir, while that of spin correlation 
remains at q—i: /2, due to the change of the relevant spin- 
orbital configuration. 

On the other hand, the Hund's rule coupling leads to 
the characteristic quantum behavior in a different way 
between spin and orbital states. In fact, the orbital 
dimer state is stabilized by the Hund's rule coupling, indi- 
cating the enhancement of quantum orbital fluctuation. 
In accordance with the robust formation of the orbital 
dimer state, we observe a development of the orbital gap, 
while the spin gap is significantly suppressed even in the 
gapped dimer state. Thus, the low-energy physics is de- 
scribed by a nearly critical spin state with the small spin 
gap on the background of the robust orbital dimer state. 
We believe that the present study is useful to explore 
novel quantum states in actual low-dimensional materi- 
als with active orbital degree of freedom. 
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